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Abstract. - In a recent letter, Martin and Zagrebnov [Europhys. Lett., 73 (2006) 1] discussed 
the thermodynamic Casimir effect for the ideal Bose gas confined in a thin film. We point out 
that their findings can be expressed in terms of previous general results for the Casimir effect 
induced by confined critical fluctuations. This highlights the links between the Casimir effect 
in the contexts of critical phenomena and Bose-Einstein condensation. 



The ideal Bose gas undergoes a phase transition in the grand canonical ensemble if the 
chemical potential /i < equals its critical value /Zc = at which Bose-Einstein condensa- 
tion takes place. This transition is accompained by density fluctuations with an increasing 
correlation length /ic) ^ {fJ-c — m)^" — \ foi' the ideal gas). 

In ref. [1] the authors consider the limit of large film thicknesses d and fixed /i of the grand 
canonical potential per unit (transverse) area fdiT, fi) of an ideal Bose gas in spatial dimension 
D ~ 3. For fi fic this implies that eventually d/£,+ ^ 1 and therefore the confining 
boundaries are subject to a vanishing Casimir force resulting from correlated fluctuations. 
On the other hand, at the critical point /x — fic, d/^+ and the fluctuations become 
long-ranged, giving rise to a Casimir force per unit area F{d,T,ii = /ic) = 2k^T/^/d^ -I- . . ., 
characterized by a universal amplitude A("'^). The cross over between these two regimes, which 
has not been investigated in ref. [1], is determined by SipdiT^jj) = ^d{T,f) — dipbuik{T, ^) — 
fsurf{T,ti) (see eqs. (2) and (9) in ref. [1]) for large d and ^xed ratio d/£_+ ^ d{—^)2. Here 
we consider the case of Dirichlet boundary conditions (BCs), although the extention to other 
cases is straightforward. According to eq. (17) in ref. [1] one has: 

^ ^ n,r— 1 n— 1 ^ ^ r— 1 

(*) E-mail: gcimbassiOmf .mpg.de 
(^)In passing we note that the amplitudes A(T) introduced in eqs. (3), (12), and (13) in ref. [1] are not 
universal, unless they are properly normalized as A = A(T)/k-BT. 
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where /3 = l/(fcBT), A = h^J j3/m is the thermal wavelength, $„(s,u) = s '^e 2n^/s^ 
and u = {-2l3fi)U/X - For d/X > 1, (A/d)^ $„((A/d)2r, u) ^ ds $„(s, u), 

with corrections decaying faster than any power of X/d, leading to [2] 

pdHMT,^i) = E EE e(u), (2) 

71—1 

so that F{d,T,n) ^ -dSipdiT, n)/dd = fcBT[2e(M) - ue'(M)]/d3. Q{u) is a continuous, neg- 
ative, and monotonicly increasing function of u providing the interpolation between the two 
cases discussed in ref. [1]: At the transition point u = one recovers eq. (13) therein, i.e., 
A = 8(0) = -C(3)/(87r). Due to e{u > 1) - e"^" for d > A and fixed ^ we recover eq. (21): 
\S(pdiT, < C)(e-\/-8/3Mrf/A)^ It is easy to verify that e{u) = Q+o,oi^) in D = 3 and with 
TV = 2, where <d%^o 

is given by eq. (6.6)(^) in ref. [3]. q is the Gaussian, (1), uni- 

versal scaling function above the bulk critical temperature, -|-, for the finite-size contribution 
l3d'^-^6Td{T) to the free energy Td{T) = dTbuik{T) + TsurfiT) + 6Td{T) per unit (trans- 
verse) area of systems the critical properties of which are captured by the 0{N) symmetric 
Landau-Ginzburg (LG) Hamiltonian [4], confined in a d x oo^~^ slab with Dirichlet BCs, 
O, O. This unnoticed connection between the results of refs. [3] and [1] holds also for periodic 
and Neumann BCs considered in ref. [1]. It is rooted in the fact that the grand canonical 
partition function for a weakly interacting Bose gas in D dimensions can be expressed, close to 
the transition point (D > 2), as a functional integral with weight e""^'*^! (see, e.g., refs. [4,5]) 
where 4'{x) is a two-component real field, S the 0(2) LG Hamiltonian 

S[cp] = J d^x |i[V</.(x)]^ -f- lr<p^x) + |[0'(x)]^| , (3) 

r = —2mfi/h'^, and g ~ ASnaK^/X^ where a is the scattering length. For the ideal gas 5 = 0, 
and S reduces to the so-called Gaussian model (defined only for r > 0), characterized by 
= r~2. 5 is a particular case of the more general 0(A^)-symmetric LG Hamiltonian 
considered in ref. [3], where the universal scaling functions 0(y+) corresponding to different 
BCs (surface universality classes) imposed on (/)(a;) have been determined analytically within 
the field-theoretical e-expansioon {e — 4 — D), as functions of the scaling variable j/+ = d/^+. 
For the ideal Bose gas j/+ ^ dr^ ~ u, as in eq. (0). Thus, the general results of ref. [3] predict 
also the Casimir amplitude Aint in a weakly interacting Bose gas (as well as for the supcrfluid 
transition, belonging to the same XY universality class N = 2). In _D = 3 and for Dirichlet 
BCs Aint — —0.022 compared to the ideal gas case A ~ —0.048. 
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{^)The correct lower integration limit is a; = 1 in accordance with eq. (6.5). 



